A Stable Pfaffian State in Bilayer Graphene 
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Here we show that the Pfaffian state proposed for the | fractional quantum Hall states in con- 
ventional two-dimensional electron systems can be readily realized in a bilayer graphene at one of 
the Landau levels. The properties and stability of the Pfaffian state at this special Landau level 
strongly depend on the magnetic field strength. The graphene system shows a transition from the 
incompressible to a compressible state with increasing magnetic field. At a finite magnetic field 
of ~ 10 Tesla, the Pfaffian state in bilayer graphene becomes more stable than its counterpart in 
conventional electron systems. 
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Ever since the discovery of the quantum Hall state at 
the Landau level filling factor v = ~, the first even- 
denominator state observed in a single-layer system, it 
has been very aptly characterized as an "enigma" 0. It 
was clear at the outset that this state must be different 
from the fractional quantum Hall effect (FQHE) in pre- 
dominantly odd-denominator filling fractions 0, U] . Un- 
derstanding this enigmatic state has been a major chal- 
lenge in all these years @ . At this half- filled first excited 
Landau level @, a novel state described by a pair wave 
function involving a Pfaffian @, 1 (or anti-Pfaffian @) 
has been the strongest candidate. More intriguing are the 
elementary charged excitations at this ground state that 
have a charge e* = e/4 and obey 'non-abelian' statistics 
[iol [TT| . Recent observation of the e* = e/4 quasiparti- 
cle charge at v = | quantum Hall state [12| has brought 
the issue to the fore [La ]. It has been suggested that 
these non-abelian quasiparticles, besides carrying the sig- 
natures of Majorana fermions [14| in this system, might 
even be useful for quantum information storage and pro- 
cessing in an intrinsically fault-tolerant manner [15j . 

Electrons in another recently discovered two- 
dimensional system, graphene Hi . display a range 
of truly remarkable behavior [17] . The dynamics of 
electrons in a single sheet of graphene, a hexagonal 
honeycombed lattice of carbon atoms is that of massless 
Dirac fermions with linear dispersion, chiral eigenstates, 
valley degeneracy, and unusual Landau levels in an 
external magnetic field [ijj ■ Theoretical studies of 
FQHE in monolayer [lj| and bilayer graphene [l9[ were 
reported earlier by us. Recent experimental observations 
of the v = | FQHE in monolayer graphene [2Cj have 
provided a glimpse of the role highly correlated electrons 
play in graphene. Given the accute interest in studying 
the properties of the v = | state in conventional two- 
dimensional electron gas (2DEG), a natural question to 
ask is how does this state manifests itself in graphene. 

For the conventional (nonrelativistic) 2DEG the in- 
compressible state at v = | has been studied numeri- 
cally for a finite number of electrons [li|. A relatively 
good (but not 100%) overlap with the Pfaffian state has 
been found. The overlap of the exact wave function of 



the finite-size systems with the Pfaffian state can be im- 
proved by varying the inter-electron potential. For exam- 
ple, by increasing the thickness of the two-dimensional 
layer 2l| , one can improve the overlap with the Pfaffian 
state and increase the excitation gap of the correspond- 
ing incompressible state. The interaction properties of a 
two-dimensional system are determined by the Haldane 
pseudopotentials [22| , which are the energies of two elec- 
trons with relative angular momentum m. The pseu- 
dopotentials at the n-th Landau level are of the form 
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where L m (x) are the Laguerre polynomials, V(q) — 
2ire / (K,£ q) is the Coulomb interaction potential in the 
momentum space, k is the dielectric constant, £ — 
(h/eB)2 is the magnetic length, and F n (q) are the form 
factors of the n-th Landau level. 

Within the framework of the Haldane pseudopotentials 
it is convenient to study the finite-size system numeri- 
cally in the spherical geometry. The size of the sphere 
and the number of single-particle states are determined 
by the parameter S, where 2S is the number of magnetic 
fluxes through the sphere in units of the flux quanta. 
The single-electron states are characterized by the angu- 
lar momentum S, and its z component S z . For the many- 
electron system the corresponding states are classifie d by 
the total angular momentum L and its z component [23| . 
For a system with iV electrons the v = ^ Pfaffian state 
is realized at 2S = 2N — 3. Here the filling factor v = | 
is defined as the filling factor of a given Landau level. In 
spherical geometry the v — \ Pfaffian state is the exact 
ground state only for a very special type of three-particle 
interaction [§[ when the three-particle interaction poten- 
tial is non-zero only if the total angular momentum of 
three particles is 3S — 3. For any two-particle interac- 
tion the v — i Pfaffian state is not an exact eigenstate, 
which makes it impossible to continuously connect the 
Pfaffian state to any exact eigenstate of the two-particle 
Hamiltonian. By varying the interaction function, i.e., 
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the pseudopotentials, the close proximity to the Pfaffian 
function with an overlap of 99% can be achieved. The 
v = i Pfaffian state is most sensitive to the lowest pseu- 
dopotentials, Vi, V3, and V5. 

For a single graphene layer the Landau level wave func- 
tions are mixtures of those for Landau levels of nonrel- 
ativistic systems; for example the first Landau level in 
graphene can be expressed in terms of zero and the first 
Landau wave functions of the nonrelativistic system [l8| . 
As a result the form-factor in a single graphene layer 
takes the form F n (q) = [L n {q 2 /2) + L n _i (g 2 / 2 )]/ 2 ' for 
n > 1 and F n=0 = L (q 2 /2) [18|. Numerical analysis of 
finite-size systems in a spherical geometry with up to 14 
electrons shows that the largest excitation gap around 
0.02 e 2 /(K£ ) occurs at the n — 2 graphene Landau level 
24] . Although the excitation gap at a finite size system 
in this case is comparable to the v = § nonrelativistic 
system the overlap of the ground state with the Pfaffian 
state is less than 0.5 at all Landau levels [13]. This fact 
shows that a single graphene layer does not have stable 
incompressible v = \ Pfaffian states. Modification of the 
interaction potential can improve the formation of the 
v = \ incompressible Pfaffian state. One such modifica- 
tion can be achieved in a bilayer graphene. 

We show here that bilayer graphene can indeed im- 
prove the stability of the v — \ Pfaffian state in graphene. 
Namely, one of the bilayer Landau levels (for a given 
valley) has a stable v = \ Pfaffian state, the proper- 
ties of which can be controlled by a magnetic field. The 
maximum overlap of the finite system ground state with 
the corresponding Pfaffian state occurs at finite values of 
the magnetic field. The v = \ incompressible state of a 
bilayer graphene is more stable than the corresponding 
state in a conventional two-dimensional system. 

We consider a bilayer graphene which consists of two 
coupled graphene layers with the Bernal stacking ar- 
rangement. Each graphene layer has two sublattices, 
say, A and B. For the Bernal stacking arrangement, the 
coupling is mainly between the atoms of sublattice A of 
the lower layer and atoms of sublattice B' of the upper 
layer. For one projection of spin, e.g., +i, the state of 
the bilayer graphene can be expressed in terms of the 
four-component spinor (tpA, ipB, V'B'i ^A') T for valley K 
and (ipB' , ipA' , ipA, V , b) T f° r valley K' . The subindices 
A, B and A', B' correspond to lower and upper layers 
respectively. The properties of bilayer graphene can be 
controlled by a bias voltage, AU, which is the poten- 
tial difference between the upper and lower layers. The 
Hamiltonian of the biased bilayer system in a perpendic- 
ular magnetic field has the form [25J 
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FIG. 1: (a) Few lowest Landau levels of a bilayer graphene, 
shown for AU = 100 meV and t = 400 meV. The two solid 
red lines belonging to different valleys, correspond to the Lan- 
dau levels where the v = \ Pfaffian state can be observed, 
(b) Overlap of the exact many-particle ground state with the 
Pfaffian function, (c) Excitation gap of the v = ^ state. The 
results are for N = 14 and 25* = 25 and the zero bias voltage. 
The black and red lines correspond to t = 400 meV and 300 
meV, respectively. Here the energy unit is e c = e 2 /k£q. 
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where t is the inter-layer hopping integral, n± — Tr x ±iTr y , 
7? = p + eA/c, p is an electron two-dimensional momen- 



FIG. 2: The excitation gap (a) and the overlap with the Pfaf- 
fian state (b) versus the bias voltage, AU. Here N = 14 and 
25 = 25, and the magnetic field is 10 Tesla. 



turn, A is the vector potential, v F s» 10 6 m/s is the fermi 
velocity, and £ = + (K valley) or — (K 1 valley). 

The discrete eigenstates of the Hamiltonian (Q} can be 
found from the following equation (25| 

"(e + i5f - 2(n + l)][(e - £<5) 2 - 2n] = (e 2 - 5 2 )t 2 , 

(2) 

where S = AU/2 and all energies are expressed in units 
of Hv F /£ . For a given value of n there are four bilayer 
Landau levels which are characterized by the index n 
and the energy e of the level. The corresponding wave 
functions can be expressed in terms of n, \n— 1|, and n+1 
conventional Landau wave functions [25] . The resulting 
form factors F nyS (q) were derived in [19| . 

The form factors and the corresponding pseudopoten- 
tials allow us to find the energy spectrum of a finite N- 
electron system in the spherical geometry [l9|, [22j • We 
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FIG. 3: The (a) excitation gap and the overlap with the Pfaf- 
fian state (b) are shown for different number of electrons: 
N = 8, 10, and 14. The magnetic field is 10 Tesla and the 
bias voltage is zero. 
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report our calculations for N = 8, 10, and 14 electron 
systems. To determine the incompressibility of the sys- 
tem we calculated the excitation gap and the overlap of 
the ground state wave function with the Pfaffian func- 
tion. We consider only one valley, for example, valley K. 
The results are similar for K'. 

For all but one bilayer Landau levels the overlap of the 
v = i ground state with the Pfaffian state is found to 
be small (< 0.6). At the same time there is one spe- 
cial Landau level (for each valley) at which the v = | 
ground state is well described by the Pfaffian function. 
This level corresponds to one of the solutions of Eq. @ 
with n = 0. At a small bias voltage, AU, the wave 
function corresponding to this Landau level is the mix- 
ture of the conventional Landau level wave functions 
with indices and 1. The wave functions of this spe- 
cial level have the form of (</>o,0,0, (t/ y/2)(f>i) , where </> n 
are n-th 'nonrelativistic' Landau functions and t is in 
units of hv F /£ . Then the corresponding form factor is 
F(q) = (L Q + (t 2 /2)L 1 )/(l + t 2 /2). At small values of 
the dimensionless hopping integral, t(£ /Kv F ), the inter- 
action within this level is similar to the one at the lowest 
Landau level of a conventional system, which does not 
show an incompressible v = \ state. At large values of 
t(£ /hv F ) the special bilayer Landau level is similar to the 
n = 1 Landau level of the conventional system and shows 
the v = i Pfaffian state. By varying the magnetic field, 
the dimensionless inter-layer hopping integral is changed 
which modify the interaction within the Landau level and 
changes the properties of the v = \ state. We present 
the numerical results only for this special bilayer Landau 
level. 

At the zero bias voltage this special Landau level has 
zero energy and is degenerate with another level, which 
has the form (0,0, 0, 4>q). In addition to this acciden- 
tal degeneracy, each level has two-fold valley degeneracy, 
which make the zero energy state four-fold degenerate. 
At a finite bias voltage this degeneracy is completely 
lifted and the special Landau level of the bilayer can be 
isolated. In Fig. 1(a) we show the lowest Landau levels 
of a bilayer at finite bias voltage. Two solid red lines 
correspond to the special Landau levels of the two val- 



FIG. 4: Trajectory of the inter-electron interaction with vary- 
ing magnetic field, shown by a solid red line in the plane 
(Vi/Vs) — (V3/V5) for the 'special' Landau level of the bilayer 
graphene. The initial point of the trajectory (at B = 0) cor- 
responds to the nonrelativistic system at n = 1 Landau level, 
while the final point (at B — 00) corresponds to the nonrela- 
tivistic system at the n — Landau level. The shaded region 
illustrates the compressible v = | state, while the blank re- 
gion corresponds to the incompressible v = \ state (Ref. [H|]). 
The crossing of the boundary between the compressible and 
incompressible states occurs at B ~ 100 Tesla for the hopping 
integral t = 400 meV. The blue dashed line shows the region 
of large overlap with the Pfaffian state (Ref. [IH). 



leys. The many-particles properties of these two levels 
are identical. Therefore, we study the v = | state for 
only one valley. 

In Fig. l(b,c) we show the magnetic field dependence 
of the overlap of the v = \ ground state with the Pfaf- 
fian state and the corresponding excitation gap. At a 
small magnetic field the dimensionless hopping integral 
is large and the system becomes similar to the conven- 
tional system at the n = 1 Landau level. With increasing 
magnetic field the properties of the system changes non- 
monotonically and the overlap with the Pfaffian state 
reaches its maximum at a magnetic field of ~ 10 Tesla 
(and for t — 400 meV). The overlap at this point is 
~ 0.92, which is a big improvement over the nonrelativis- 
tic system (~ 0.75). The dimensionless hopping integral 
at this point is t(£ /Hv F ) ~ 4.89. 

At a large magnetic field the system is close to the 
n = nonrelativistic Landau level, the overlap with the 
Pfaffian state becomes small and the v — | state is fi- 
nally compressible. This dependence on the magnetic 
field opens up interesting possibilities to investigate the 
stability and appearance and disappearance of the v = h 
Pfaffian state in a single bilayer Landau level. Although 
the Pfaffian state becomes unstable only at large mag- 
netic fields, this property strongly depends on the value 
of the hopping integral. At smaller hopping integrals 
the magnetic field range of stability of the Pfaffian state 
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shrinks. For example, at t = 300 meV the Pfaffian state 
is expected to be unstable at B ~ 40 Tesla (see Fig. 1). 
Another parameter which controls the properties of the 
graphene bilayer is the bias voltage. Although the bias 
voltage modifies the bilayer wave functions, we found 
that the overlap of the ground state with the Pfaffian 
state and the excitation gap have weak dependence on 
the bias voltage within a broad range of AU (see Fig. 2). 
The overlap monotonically decreases with increasing AU, 
which suppresses the overlap by only a few percent. The 
large excitation gap and the large overlap observed for 
different system sizes are shown in Fig. 3. 

The stability and the strength of the Pfaffian state can 
be also analyzed in terms of the general dependence of 
the pseudopotentials, V m , on the relative angular mo- 
mentum, m. We characterize the interaction potential 
of the bilayer graphene by two parameters: V4/V5 and 
V3/V5 13]. These parameters depend on the magnetic 
field. By varying the magnetic field, this dependence can 
be shown as a line in the (V1/V5) — (V3/V5) plane (Fig. 4). 
That line connects the initial point at B = to the fi- 
nal point, corresponding to large magnetic field, B = 00. 
The v — i bilayer graphene system at the initial and 
final points are identical to the conventional systems at 



the first (n = 1) and zero (n — 0) Landau levels, respec- 
tively. In Ref. |13J the region of the compressible v = i 
state and the region of strong overlap with the Pfaffian 
state were identified (see Fig. 4). With increasing mag- 
netic field, the v = h bilayer graphene system transforms 
from a v = | nonrelativistic state (at small values of B) 
to a more stable incompressible state with large overlap, 
and finally to a compressible state (at a large magnetic 
field). For the hopping integral t — 400 meV, the tran- 
sition from the incompressible to a compressible v = h 
state occurs at B ~ 100 Tesla. 

In conclusion, a stable incompressible v = i Pfaffian 
state can in fact be observed in a bilayer graphene only 
at one Landau level. The properties of this state strongly 
depend on the value of the magnetic field. With an in- 
creasing magnetic field, the v = \ state transforms from 
an incompressible state at a small magnetic field to a 
compressible state at a large magnetic field. At inter- 
mediate values of the magnetic field, B ~ 10 Tesla, the 
v = \ state becomes more stable than the corresponding 
state in a conventional two-dimensional electron system. 
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